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Abstract 

The phase diagram of the unconventional superconductor UPts is explained 

under the long-standing hypothesis that the pair wavefunction belongs to the 

Eig representation of the point group. The main objection to this theory 

has been that it disagrees with the experimental phase diagram when a field 

is applied along the c-axis. By a careful analysis of the free energy this 

objection is shown to be incorrect. This singlet theory also explains the 

unusual anisotropy in the upper critical field curves, often thought to indicate 

a triplet pair function. 
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Unconventional superconductivity is a state of matter under intense discussion at the 
present time, in both high-Tc materials and the older heavy fermion superconductors. In 
this latter class of materials the most studied and best characterized is UPts. The d-wave 
Eig state was originally proposed as the pairing symmetry on microscopic grounds and 
has had good success in explaining a large class of experiments. It posits a two-component 
gap function (in contrast to the d-wave states believed to be relevant to high-Tc) which 
transforms as {k^kz, kykz) with corresponding line nodes where the Fermi surface intersects 
the plane kz = and point nodes where it intersects the line k^ = ky. Evidence for this 
specific pattern of nodes comes from ultrasound and heat conduction experiments. 
Eig also explains the pressure dependence of the phase diagram. Eig, along with other 
two-dimensional representations of Dq^ has a two-component order parameter (OP). This 
leads to a number of unusual predictions which have been confirmed by experiment, for 
example the split transition in specific heat measurements and the kink in the lower 
critical field curve A two-component OP is usually (though not always ^j) accepted for 

UPt3. 

In spite of the fact that has the proper nodal structure and number of components, a 
number of alternatives have been proposed. The objection usually given is that Eig cannot 
explain the observed phase diagram in the field-temperature (H-T) plane when H is along 
the c-axis 0. A second objection to Eig, a singlet theory, is that the upper critical field 
curve Hc2x{T) for H in the basal plane crosses the curve Hc2z{T) for H perpendicular to 
the basal plane |jlO| and that this is characteristic of triplet theories. |[ll|] In this letter we 



show that both objections are unfounded. 

It has been clear for several years that the Eig theory predicts correctly the exceedingly 
unusual phase diagram in the H-T plane when H is in the basal plane. ||12| There are 
three superconducting phases meeting the normal phase at a tetracritical point. Two of 
these, the A and C phases, are conventional distorted Abrikosov lattices formed by one of 
the two components of the OP. The third, the B phase, consists of two interpenetrating 
lattices, one formed by each component. The phases are separated by second order phase 



boundaries whose properties (such as the specific heat jump ACy) may be calculated. These 
conclusions and the conclusions of the present paper follow from the free energy density for 
the Eig theory: 

/ = ao{T - TJIr^^f + ao{T - T,)|r/,|2 + /5i(r7 ■ rj*f + P^lv ■ vl' (1) 

i,j=x,y i=x,y 

+ {aoe AT) {he/ 2e) ^ (|Ar/.f - lAr?,!') 

i=x,y 

+ a.Hlr] ■ ri* + a^{H^ + H^)ri ■ ij* + ad\H ■ ■q\^ . 

Here 77 = {rjx.rjy) is the two-component order parameter, and Ki, K2, K3, K4, ao, Pi, P2, 
a^i CLz) cid and e are constants and AT = Tj. — Ty. The D's are momentum operators: = 
—id/dx + {2e/hc)Ax and similarly for Dy and D^. Here A is the vector potential and — e is 
the charge on an electron. The coupling of the staggered magnetization to rj is responsible 
for the temperature splitting AT. The existence and need for the term proportional to e, 
which represents the coupling of the supercurrent to the staggered magnetization, and the 
terms proportional to was first stressed in the context of a threee-component model. |TB 



To obtain Hf.2 we need only consider the terms quadratic in 77 in Eq. (0). We minimize 
these terms following the Euler-Lagrange perscription of 5F = where F = J d?xf . When 
the field is in the basal plane this procedure decouples into one d.e. for rj^ and one d.e. for 
rjy. Hence we obtain two separate equations for Hc2- For appropriate values of the constants 
these two curves will cross creating the well-known kink in the upper critical field curve. 
Hence the A and C phases correspond to ry ~ (1, 0) and rj ~ (0, 1) in our theory. The d.e.'s 
for rjx and rjy in this case both have the same form as the corresponding equation in the more 
familiar single component Ginzburg-Landau theory. Hence in either phase the component 
of rj that is non-zero will form an Abrikosov lattice of the usual kind. The lattice will be 
distorted however because of anisotropy in the gradient terms due to the inequality of the 
K^s. In the B phase we have two flux lattices: one formed by rjx and the other formed 
by rjy. The zeroes of these flux lattices need not coincide, however. Hence, we must 
introduce the extra degree of freedom of the offset vector of the two flux lattices, and we 



must minimize the free energy with respect to this variable. We find that the free energy is 
minimized for an offset vector which is one-half of a flux lattice basis vector. |]I4| 

When H is in the z-direction, the problem of minimizing the free energy is far more 
difficult to solve. The eigenfunctions of the linear problem, which we shall call (f)nk, can 
be found numerically but they are complicated. Because the linear Hc2 equations do not 
separate into separate equations for rjx and rjy, the 4>nk have both x and y components. 
When the OP expanded in terms of the eigenfunctions: t] = J2nkCnk4>nk, the free energy 
-F is a polynomial in the coefficients: F = F{cnk)- Here n is a level index (no longer 
a Landau level index) and k is the momentum in the y-direction. Part of the argument 
against Eig runs as follows. At Hc2, some of the Cofc become nonzero. The fourth order 
terms in F are complicated. If we take, for example, r] = cok4>ok + Cik4>ik, then terms of 
the form | Cq^ P c^k^ik appear. The c^k produce, effectively, a linear term in the cik- It is 
then concluded that no second transition occurs below Hc2 in this theory, in conflict with 
experiment. 

This argument is, however, not correct. The actual Hilbert space for the functional F is 
infinite-dimensional and a careful analysis of all the possibilities must be carried out. The 
energy levels labelled by the integer n > are highly degenerate, the eigenvalues being 
independent of which may take the value of any allowed wavevector k = 271 ^ integer/L^, 
where Ly is the length of the sample in the y-direction. The energy of the OP configuration 
represented by (t)nk is independent of k. The minimization of F leads to some of the Cofc 
becoming nonzero at ifc2 with the formation of the usual hexagonal lattice: Cnk ~ '^n.o(-f^c2 ~ 
Hy/'^Ck- Let 271 / q be the periodicity of the fiux lattice in the y direction. Then Cfc=0 unless 
k = mq, where m is an integer. As usual Ck = for m = even(odd). A dangerous fourth 
order term in F has the form: l3klk2k3ki^oki^ok2^ok3^ik4- Momentum conservation implies that 
the coefficient Pklk2k3k4 is only nonzero if fci + ^2 + ^3 + ^4 = 0. For an interpenetrating 
lattice where the offset vector is one-half of a flux lattice basis vector, ki, k2, and k^ are 
integer multiples of q, whereas k^ is half an odd integer times q. Thus the fc's never sum 
to zero and all dangerous terms vanish. The cik for the second lattice never appear in flrst 
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order or, by a similar argument, in third order. The second lattice appears by a second order 
transition in the Eig theory for all directions of the applied field. This is in agreement with 
experiment and in conflict with previous theoretical conventional wisdom. The transition 
breaks the flux lattice symmetry because the lattice now has a basis. 

We have plotted the phase boundaries obtained by minimizing the free energy in the 
following approximation. The eigenvalues of the linear Hc2 operator are obtained by a 
truncation of the inflnite matrix. The lowest eigenvalue, which is a function of H and T, 
gives the Hc2 curve. The next lowest eigenvalue gives a bare inner transition line. This 
must be corrected by an effective fleld term because the existing lattice lowers the transition 
temperature of the new one. This correction involves only one coupling constant which is 
obtained by fltting to the data. The result is shown in Fig. 0. We have not attempted to 
flt the data for T < OAT^ since the linear temperature dependence of the flrst two terms in 
Eq. (|l]) breaks down there. 

There is no tetracritical point for H = Hz; this is due to level repulsion. We regard 
this as a virtue of the theory, because the experimental data show that to call the phase 
diagram isotropic is an exaggeration. The Hc2 curve for H = Hz does not have a kink, 
only a flat region well reproduced by the theory, and the data are consistent with only two 
superconducting phases, as the present theory predicts for this fleld direction. Another part 
of the argument against Eig has been that flne tuning of parameters is required to flt the 
data. Our flt does not does not require any flne-tuning. We flnd that the phase diagram 
for both fleld directions can be flt by the same set of parameters, and the only numerical 
coincidence which arises when this is done is that K2 ~ K^. This is actually a consequence 
of approximate particle-hole symmetry and the fact that it comes out of the flt is further 
evidence that the overall picture is correct. 

To understand the directional dependence of Hc2, it is flrst necessary to discuss the 
magnetic susceptibility of UPts. This issue is complicated by the fact that all renormaliza- 
tions involved are not well understood. Since UPts is a Fermi liquid, however, the starting 
point must be the single-particle states calculated in band theory, which accounts very 
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well for the Fermi surface. |jT5|] The states near the Fermi surface are predominantly de- 
rived from uranium 5f orbitals with j = 5/2. In the isolated atom, these would be 6-fold 
degenerate. In the hexagonal crystal field, there is an effective Hamiltonian of the form 
Hcrystaifieid = Bh^jl— ]{] + 1)/3), whcrc Bh is a constant. This splits the six-fold degenerate 
state into three doublets at the F point: jz = ±5/2, = ±3/2, and jz = ±1/2. There is 
also an even-odd splitting from the fact that there are two U atoms per unit cell. The six 
bands constructed from these states cross the Fermi energy, and the crystal field splitting 
is of the same order as the bandwidth. The average occupation of the 5f level is between 2 
and 3. If we apply a magnetic field, there will be both a Pauli (intraband) and a Van Vleck 
(interband) contribution to the susceptibility. The former is of order {gefffiB)'^N{eF), while 
the latter is of order {gegUBY / \ Bh\- Here g^ff is an effective g-factor for the coupling of the 
field to the total angular momentum of the band or bands involved. It is a dimensionless 
number of order unity. The Lande factor for £=3, s=l/2 , and j=5/2 is 6/7. 
The Van Vleck susceptibility is given by 

X.. = 2n,|E^^^^^i§^/.(l-/.). (2) 

Here /«, //3, Ep are occupation factors and energies of the states a and (3. In view of 
the greater multiplicity of the interband transitions, we expect the Van Vleck susceptibility 
to be very important - indeed it may dominate the total. If H is along the c-axis, then the 
relevant matrix element (with ^ = 1) is: 

I < a\Lz + 2Sz\l3 > p = {m/m)3%,p. (3) 

In the approximation that states of different jz do not mix (negligible intersite interactions), 
then the perturbation introduced by H is diagonal, and the occupation factors then imply 
that the Van Vleck susceptibility is zero for this direction. If H is in the x-direction, the 
corresponding expression for the square of the matrix element is 

I < + 2S,\I3 > p = (36/49) (5/2 - j,)(5/2 + + 1) (4) 
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if the states a and ^ differ by one unit of jz and is zero otherwise. The Van Vleck suscep- 
tibihty comes from four distinct pairs of states:(j^=-5/2,-3/2), (-3/2,-1/2), (1/2,3/2) and 
(3/2,5/2), whenever one of the pair is occupied and the other unoccupied. The Pauh con- 
tribution to Xxx^ on the other hand, comes only from the pair (-1/2,1/2) when this state is 
occupied. 

Summing up these considerations, we expect that Xzz will be dominated by the Pauli 
contribution, while Xxx will be dominated by the Van Vleck contribution. There are two 
effects which can vitiate these conclusions. Intersite effects will mix states of different j^, 
and this will modify this ionic picture. The interaction effects give rise to the large Fermi 
liquid enhancement of the susceptibility, which comes chiefly from the mass term. This is 
expected to affect Pauli and Van Vleck terms alike. [|1^] Experiment confirms these theo- 



retical arguments. It is found that Xxx is considerably larger than Xzz at all temperatures, 
in accord with the expectation that the Van Vleck contribution is large. The temperature 
dependence of Xxx(^) is anomalous, with a peak at T=15 K. [|T^ This peak is absent in the 
smooth curve for XzzK^\ and in the the specific heat Cy{T\ |[TB| This is consistent with 
the idea that the physical origins of Xzz and Xxx are different, and that the density of states 
at the Fermi level determines Xzz but not Xxx- Thus experiment confirms the theoretical 
picture. 



The importance of these considerations for the superconducting state is simple. |T9 
Superconductivity affects the Pauli susceptibility in a drastic fashion. For a singlet state 
such as Exg-, the Pauli term xfjiT) is reduced to zero at zero temperature because it takes 
a finite amount of energy to break a pair and magnetize the system. Superconductivity 
should have no effect at all on the Van Vleck term. The difference in free energies between 
the normal and superconducting states in a field is Fmagnetic = ^XijHiHj. Here 

^Xij = Xij ~ Xij where xfj and Xij are the Pauli susceptibilities in the superconducting 
and normal states, respectively. Hence we expect a field along the c-axis to have the largest 
effect on superconductivity. Near Tc, the slope of Hc2, larger in magnitude for H in the 
z-direction, is determined by the terms in F which are linear in H. The different slopes 



reflect the anisotropic coherence length and are not directly related to the susceptibility. 
As H increases, the terms become more important and cause Hc2{T) to curve down. 
The anisotropy in the Pauli susceptibility then causes Hc2z to curve more strongly with the 
result that the two curves cross. To implement this quantitatively, we note that the change 
in the susceptibility is quadratic in rj near T^. The expression for Fmagnetic which results is 
precisely the last three terms, proportional to H^, in Eq. (|I|). The resulting fit is shown in 
Fig. |. 

What these arguments show is that the peculiar anisotropy of the upper critical field 
is evidence for a singlet superconducting state, such as Eig. This is in sharp contrast to 
previous arguments that the anisotropy points to a triplet state. These arguments were 
based on the idea that the observed anisotropy in the total susceptibility is also reflected in 
the Pauli term, that is Xxx ~ '^Xzz- According to the above arguments, this appears unlikely. 

We conclude that the Eig theory can account for two crucial aspects of the phase diagram 
of UPts: the existence and shape of the inner transition line for H along the c-axis, and 
the peculiar anisotropy of the upper critical field. This removes the major objections to 
this theory, which otherwise gives a good account of the low temperature thermodynamics, 
including the position of the gap nodes, the tetracritical point, and the dependence of the 
phase boundary positions on applied pressure. 

We are very grateful to D. L. Cox for useful discussions and acknowledge the support of 
the National Science Foundation through grant no. 9214739. 
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FIGURES 

FIG. 1. Phase diagram when the field is in the z-direction. The hnes are the theoretical fits to 
Hc2 (solid line) and the inner transition (dashed line). The data points are from ultrasonic velocity 
measurements and are taken from Ref. Fig. 3. 

FIG. 2. Shows the crossing of the Hc2 line when the field is the basal plane (solid line) and 
the Hc2 line when the field is in the z-direction (dashed line). The data points are from ultrasonic 
velocity measurements and are taken from Ref. Fig. 3. The diamonds are for the case when 
the field in the basal plane {H \\ ab) and the crosses are for the case when the field is in the 
z-direction {H || c). 
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